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Abstract 

We discuss the rapidity distribution of produced jets in heavy-ion 
colhsions at LHC. The process ahows one to determine to a good 
accuracy the value of the impact parameter of the nuclear collision 
in each single inelastic event. The knowledge of the geometry is a 
powerful tool for a detailed analysis of the process, making it possible 
to test the various different elements which, in accordance with present 
theoretical ideas, take part to the production mechanism. 
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1 Introduction 



Within a few years it will be possible to observe at LHC the first events with 
heavy ions colliding at energies of several TeV in the nucleon-nucleon cm. 
system. The interest in the process is justified by the expectation that a new 
phase of matter, namely quark-gluon plasma, could appear as a result of the 
large energies - high densities which will be reached accelerating heavy ions at 
LHC energies On the other hand the large number of fragments, resulting 
from such large structures, as say Pb nuclei, interacting with so large energy 
in the cm. system, makes the analysis of the events rather difficult and 
the search of unambiguous signatures puzzling. The identification of a few 
features in the event which can be easily understood and whose dynamics 
is, on the contrary, rather transparent would be in this respect an important 
handle. The production of jets is, in our opinion, a process which has good 
possibilities to be understood in a rather satisfactory way, also at such an 
extreme regime as in the case of heavy ion collisions at LHC [0. In fact, in 
the present paper, we discuss the mechanism of production of large kt jets and 
of minijets in heavy ion collisions at LHC energies. The main observation is 
that one has the possibility of disentangling the geometrical features and to 
obtain as a consequence distributions with universal characteristics. When 
geometry is decoupled one has in fact the possibility of constructing various 
distributions of jets, produced at different rapidities, with different cut-off 
values etc., which are related one to another in a simple way, that is easily 
worked out in perturbation theory. This universality property, on the other 
hand, represents a regularity feature which is easily tested experimentally and 
which is rather independent on the details of the actual theoretical approach. 



2 Discussion 

The inclusive cross section to produce jets at large kf is expressed, in the 
QCD-parton model, by the factorized form 

^^■^ G{x^)G{xj,)^ (1) 



dxAdxB(Pkt d?kt 

where G{x) is the parton distribution, depending on the fractional momen- 
tum X and on the scale /ct, and is the partonic cross section. By intro- 
ducing a lower cutoff in the transverse momentum of the produced jets k^ 
one defines the integrated cross section: 

= I -J — dxAdxsd^kt (2) 
J sxAXB>Hk^Y dxAdxsd^kt 
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The lower limit for the integration in xa and xb is obtained by the condition 

sxaXb > 4(A;^)^. The differential cross section as a function of xa (or as a 
function of xb) results from the corresponding integration on xb (or on xa)- 
For Xa one has to integrate xb with the limits 1 > xb > '^{kf)'^ / sxa- The 
expression of the differential cross section as a function of xa is therefore 
written as 



The relation corresponds to a flux of projectile partons, at a given value of 
momentum fraction xa-, interacting with the target. The integral represents 
the target scattering centers which are seen by a projectile parton with mo- 
mentum fraction xa- One may notice, by looking at the integration limits of 
the integral on xb-, that when xa is large the number of scattering centers 
is also large, while at small xa the number of scattering centres drops to 
zero. The size of each scattering centre is set by the value of the partonic 
cross section, and, although it depends on the value of the cut-off k^, for 
sensible values of /c^ , i.e. for values of the order of a few GeV , the size of the 
scattering centres is very small as compared to the hadronic dimension. 

When considering large nuclei interacting at LHC energy the number of 
scattering centres which is effective for producing jets with momentum trans- 
fer larger than a few GeV is a rather large number. The scattering centres 
are distributed in the volume occupied by the target nucleus. One may figure 
out the interaction process by integrating the hadronic parton distributions 
with the limits discussed above and by multiplying the result by the atomic 
mass number. The resulting quantity is the number of scattering centres as 
a function of the momentum fraction of the projectile xa- One may then dis- 
tribute the scattering centres uniformly in the volume of the target nucleus 
and take the projection in the transverse plane. In correspondence of each 
scattering centre, projected in the transverse plane, one can draw a dot whose 
size is the parton-parton interaction cross section integrated with kt > k^. 
The result is shown in fig.l for collisions with 7TeV in the nucleon-nucleon 
cm. system, 5GeV as a value of /cf and Pb as a nuclear target. The two 
figures show the distribution of scattering centres for two different values of 
Xa, Xa = 10^''' and xa = 10^^, corresponding to the case of mini-jets pro- 
duced in the central and in the forward rapidity region. The scale in fig.l is 
the fm and the size of the dot representing each scattering centre has been 
obtained by computing the partonic cross section in Born approximation. In 
the actual forward rapidity case the number of possible scattering centres 
which are effective is of the order of 20.000. 

Fig.l points out in a spectacular way one of the most relevant problems 




da. 



da. 



(3) 
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one has to deal with when considering production of jets in high energy nu- 
clear collisions. Namely the large number of elementary interactions which 
take place in a single inelastic event and the geometrical features which result 
from the localization of the elementary interactions. The first observation is 
that when the number of scattering centres is very large, there may be over- 
laps in a few cases. As a consequence the effective number of scattering 
centres is smaller with respect to the number that has been estimated simply 
by multiplying the parton distributions with the atomic mass number. To 
keep this feature into account one should consider the gluon-fusion process 
[0], which reduces the parton population accordingly. One can then observe 
that the most peculiar feature of the overall interaction is represented by the 
large number of elementary partonic collisions, which take place at differ- 
ent points in the transverse plane, inside the overlap region of the matter 
distribution of the two interacting nuclei. One may in fact expect that the 
parton population of the projectile is also distributed in transverse space in 
a way similar to the distribution of the scattering centres shown in fig.l. The 
hard, or semi-hard, nuclear collision is then obtained by superposing two 
pictures similar to those shown in fig.l, one picture representing the distri- 
bution of projectile partons and one representing the distribution of target 
partons. Both projectile and target partons are localized in transverse plane 
in a region whose size is set by the cut-off k1 and which is of the order of 
the size of the partonic interaction cross section, as in the case of the dots in 
fig.l. Projectile and target partons interact when the regions where they are 
localized in transverse plane overlap. The interactions are then localized in 
the overlap region of the matter distribution of the two nuclei and the whole 
process depends strongly on the nuclear impact parameter as a consequence. 

There are two qualitatively different kinds of multiple partonic interac- 
tions: 

• different pairs of partons interact at different points in transverse plane, 
that one may call disconnected interactions [Q, and 

• a given projectile parton may overlap at the same time with two or 
more target partons, that corresponds to a rescatterings ^ process. In 
this respect one may notice that there is no need to be in the regime 
where parton fusion is an important effect in order to have a sizeable 
rescattering probability. 

The final state produced by a single elementary partonic interaction is rep- 
resented by two or more scattered partons. The simplest possibility is to 
have just a Born level process in such a way that one has a one to one corre- 
spondence between initial state partons and jets observed in the final state. 
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Further jets produced in the single hard interaction are then to be regarded 
as corrections of order to the lowest order process. 

The quantities which are needed for the actual description of the nuclear 
hard interaction process are the (non-perturbative) initial state partonic nu- 
clear distribution and the (perturbative) interaction probability. Because of 
the multiple partonic interactions, the information which is needed on the 
initial state is however much more detailed with respect to the nuclear par- 
ton distributions, usually considered in inclusive processes. As one realizes 
looking at fig.l, one needs in fact to know not only the number of partons 
with given momentum fraction but also their position in transverse plane. 
In addition, given the number of partons with given momentum fraction 
and at a given position in transverse plane, one needs to know the parton 
population at a different transverse point and with a different x value. The 
non-perturbative input to the process is in fact represented by the whole 
nuclear multi-parton distribution , which corresponds to a much more de- 
tailed knowledge on the nuclear parton structure than presently available. 
One needs in fact to know the whole set of all multi-parton correlations, 
which are independent quantities with respect to the parton distributions 
usually considered. 

It seems however reasonable that when dealing with many partons local- 
ized in different space-time regions, as in the case of the parton population 
of a large nucleus, correlations are not a major effect. When correlations 
are completely neglected the parton population is described by a Poisson 
distribution. The probability to have a configuration with n partons with 
coordinates xi, bi, . . . x„, b„ is therefore given by 

r(xi,bi,...x„,b„) = ^(^i^bO-. /^(a:.,b„) ^_j^(.,b)...^, 

n! 

In this simplest case the whole distribution is expressed by means of the 
average number of partons with coordinates x, b, namely D{x, b). The con- 
nection with the usual parton distributions is / D{x,h)(Pb = G{x), where 
G{x) is the nuclear parton distribution usually considered in large kf pro- 
cesses. The whole multi-parton distribution is singular in the infrared region 
and the x— integration in the normalization factor is done by keeping into 
account the cutoff k^. 

Given a configuration with n projectile partons and m target partons, 
the interaction probability is obtained from the elementary probability of 
interaction pij, where i and j label a given projectile and a given target 
parton respectively. Then 1 — pij is the probability of no interaction for the 
pair i,j and Hi^j'^ll ~Pij) represents the probability of no interaction for the 
configuration with n projectile partons and m target partons. The quantity 



5 



which one needs for the cross section is the probabihty to have at least one 
interaction 0, namely P(n, m) = 1 — Y[i 'P{l ~Pij), which is in fact equal to 
the sum over all interaction probabilities. 

The hard cross section aH is obtained by integrating on the impact param- 
eter (3 the probability to have at least one hard interaction. The probability 
to have at least one hard interaction in the nuclear collision is the result of 
the sum over all partonic configurations of the projectile and of the target nu- 
clei, multiplied by the probability to have at least one hard interaction in the 
given partonic configuration. The hard cross section is therefore expressed 
as i 

(Th= f d^PY.^n-Tm-P{n,m) (5) 

n,m 

The relation above is too complicated to be read off and computed di- 
rectly, one can obtain a closed analytic form if one takes the drastic attitude 
of neglecting all parton rescatterings. One may however derive a few simple 
relations out of Eq.(5). If one works out the average number of parton colli- 
sions (u) ^ the result is the single scattering expression of the perturbative 
QCD parton model: 

(z/) X an = [ G{xA)G{xB)4^dxAdxb(fkt (6) 

JsxAXB>4:{k^f d fit 

The relation shows that a rather general non-trivial feature, namely the AGK 
cancellation, which is expected to hold as a consequence of the AGK cutting 
rules when multiple interactions are present, is implemented explicitly in 
the actual case. 

A quantity of interest is the average number of wounded partons [§], 
namely partons which have undergone at least one hard interaction. Those 
are in fact the initial state partons that are observed as jets in the final state. 
In the simplest case, when the elementary parton-parton interaction is repre- 
sented as a two to two process, all jets observed in the final state are wounded 
partons. The average number of wounded partons of the projectile nucleus 
A, with a given momentum fraction xa in a nuclear collision with impact 
parameter /5, W{xa), can be derived from Eq.(5). The resulting expression 
is: 

W{xa) = J d'^bD{xA,i> — (3) ■ — exp[— J D{xB,h)a{xA, XB)dxB]^ (7) 

a{xA, Xb) is the partonic cross section, integrated on the momentum transfer 
kt with the cut off and xb is integrated within the limits 4(fc^)^/sxA < 
xb < 1- 

One may notice that, while the initial state parton distribution is singular for 
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small X and k^ values, the average number of wounded partons is on the con- 
trary a smooth function in the infrared region. In fact the expression above 
results from the product of the average number of A-partons, D{xa-, b — 
and the interaction probability, namely | l—exp[— J D{xb, h)a(xA, XB)dxB] \- 



At small xa the number of projectile partons is very large, in that limit the 
interaction probability however goes to zero. The interaction probability in 
addition has a finite limit, actually the black disk limit, when — > 1 or 

when /c^ ^ 0. 

A further observation is that, when the interaction probability is small, if the 
dependence of the average number of partons on x and on b is factorized, 
the expression for the average number of wounded partons is also factorized 
in the same variables. More explicitly, if one writes 



where T(b) is the nuclear thickness function normalized to one, one obtains, 
in the limit of small interaction probability. 



W{xa) = G{xa) / G{xbMxa, XB)dXB X / d%TA{h - l3)TB{h) (9) 



The average number of wounded partons W{xa)-i in a nuclear collision with 
impact parameter /9, is therefore proportional in this case to the overlap of 
the matter distribution of the two interacting nuclei / d^6r^(b — /5)rB(b). 

Given the explicit form of the cross section, one may work out all higher 
moments of the distribution in the number of wounded partons. The re- 
sulting expression is more and more involved when higher and higher mo- 
ments are considered. An important simplification is however obtained, in 
the distribution of the wounded partons of the projectile, if one makes the 
hypothesis of neglecting all interactions where different projectile partons hit 
the same target parton. The set of events of this kind is rather small when 
compared with the whole set of possible events. When this simplification is 
made the distribution in the number of wounded partons, at a given value of 

turns out to be a Poissonian and it is therefore determined completely by 
its average value. A characteristic feature of the interaction is therefore that 
the fluctuation of the distribution in the number of wounded partons is less 
and less important, as compared with the average value, when the average 
value grows. It is precisely this feature that allows one to disentangle the 
geometrical aspect of the interaction event by event. 

When the distribution of the wounded partons is known one can easily 
simulate the production of jets in a collision involving heavy nuclei at LHC 
energy. In flg.2 one may see how the number of produced jets is distributed 




(8) 
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in rapidity in a central Pb — Pb collision with 7TeV energy in the nucleon- 
nucleon cm. system. In fig. 2a one shows the distribution of jets produced 
with a transverse momentum larger than 20GeV, in fig. 26 the cutoff is 
lowered to lOGeV and in fig. 2c = 5GeV. In the figures the histograms 
represent the simulated event. The continuous line is the average number of 
wounded partons, as it results after summing projectile and target wounded 
partons, as expressed in Eq.(7) using rapidity as a variable. The dashed line 
is the same average number when using rather the approximate expression 
in Eq.(9), namely neglecting semi-hard parton rescatterings. 
One may notice that while the average number grows the fluctuation around 
the average grows much less rapidly. The effect of the unitarization of the 
hard parton-nucleus interaction, namely the difference between the contin- 
uous and the dashed line, is also shown to be a rather sizable effect when 
k^ = 5GeV. 

Since the fluctuations are relatively less important when the average num- 
bers are large, one may try to determine the value of the impact parameter 
by counting, in a given event, all jets produced above some lower cutoff k^. In 
fig.3 we show the distribution in the total number of jets produced in central 
collision with different cutoff values {20GeV fig. 3a, lOGeV fig.36 and 7GeV 
fig. 3c). 

In the case of a 7GeV cutoff the fluctuation in the number of produced 
minijets is 2.5%. The overall number of produced minijets in a central colli- 
sion with a 7Gev cutoff is however larger than 1500, and the actual possibility 
of detecting such a large number of minijets is out of reach. A quantity which 
is hkely to be more accessible experimentally could be the overall energy car- 
ried by all the produced minijets. Same numbers of produced minijets can 
carry different amounts of energy. The relative fluctuation in energy is there- 
fore larger with respect to the relative fluctuation in the number of produced 
minijets. Anyhow with a 7GeV cutoff the fluctuation in energy turns out to 
be less than 4%. If one lowers the cutoff to 5GeV the fluctuation in the en- 
ergy carried by all the produced minijets is close to 2.5%. A similar variation 
in the average energy of produced jets is obtained by varying the value of the 
impact parameter by an amount sizeably smaller than Ifm. It seems there- 
fore plausible that the measure of the total energy carried by jets, produced 
with relatively small momentum transfer, will allow the determination of the 
value of the impact parameter of the nuclear collision with an uncertainity 
smaller than 1 fm. 

A possible procedure in classifying events with jets in heavy ion collisions 
could be therefore the following: 

In each single event one may measure the total energy of the produced jets 
as a function of kf. Since the distribution is narrower and narrower when k9 
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is smaller and smaller, the energy carried by the jets with relatively small 
allows one to estimate the value of the impact parameter of the collision. 
One may then collect all events with the same impact parameter f3. At each 
rapidity value and for a given cutoff one has therefore a distribution in 
the number of observed minijets. One expects that all these distributions 
are universal, to a large extent. One expects in fact that the distributions 
are the same at different rapidity values and with different cutoffs, the only 
difference being in the average values. The dependence of the average values, 
as a function of rapidity and cutoff, is obtained through pQCD. 

As an example one may write the expression for the average number of 
produced minijets, at a given rapidity y, for all events with the same impact 
parameter f3. When is not too small, say it is larger than 10 — 15GeV, 
semi-hard rescatterings can be neglected and one obtains: 



^ = I d'bTAih - p)TB{h) x[gA{y)Iy,^ygB{y'My-y',ktW 



dN 
dy 

+9B{y) ly^y' QAiv'W - y, my'] (10) 



here g{y) = xG{x) and rapidities and fractional momenta are connected by 
the usual 2 — > 2 kinematics. 

The dependence on y and on the cutoff k^, which is described by the term 
in square parenthesis, is obtained by the usual pQCD parton model inputs. 
Geometry is factorized and it enters in the first term in Eq.(lO) only. All 
different averages of distributions with the same impact parameter are related 
one to another according with pQCD as shown in Eq.(lO). When the impact 
parameter is changed the distributions are obtained, from the corresponding 
ones measured previously, simply by rescaling the averages with the value of 
the overlap function. 

One expects in addition that if one plots the relative rates of events with 
different numbers of minijets, with kt larger that some fixed value k^ and at 
the same rapidity, one obtains a quantity which tests the overlap function, 
namely which is sensitive to the actual distribution in space of the nuclear 
partonic matter. As an example in fig. 4a we plot the relative rates of events 
with jets at 7/ = and k^ = 20GeV in the case of a Woods-Saxon distribution 
in space of the nuclear partonic matter (continuous histogram) and in the 
case of an hard-sphere distribution (dottet histogram). The difference is 
much more evident if one plots the ratio of the two histograms (fig.46). 



9 



3 Conclusions 



A large number of jets are expected in the case of nuclear collisions at LHC 
The global features of the events are understood in terms of the pQCD parton 
model, once the geometrical features of the interaction are taken into account. 
The picture of the interaction just described is easily tested: One may first 
look at central collision events only, by taking a veto trigger in the forward 
direction. All central events are expected to show a universal behavior in the 
distributions of jets produced at different rapidity values and with different 
cutoff values. The jets observed are in fact the wounded partons, namely 
the initial state partons which have undcgone at least one hard interaction. 
The expectation is that the different distributions in the number of wounded 
partons, and therefore in the number of observed jets, are close to a Poisson 
distribution when the impact parameter of the nuclear collision is fixed. At 
fixed impact parameter all different distributions in the number of observed 
jets are therefore equal. The only difference is in the average value, which is 
calculation in the pQCD parton model in a staightforward way. 

One can then proceed estimating the impact parameter event by event. 
That can be done by looking at the total energy carried by the jets produced 
in each single event at small k^. The value of the energy, compared with the 
corresponding value in the case of a central collisions, gives the estimate of 
the size of overlap of the two interacting nuclei. Within the simple approach 
to the nuclear interaction actually described, collecting all events with the 
same overlap one obtains distributions which have to be the same as those 
observed in the case of central coUisions, after rescahng the averages with 
the corresponding value of the overlap function. 

The simple picture of the interaction just described is easily tested ex- 
perimentally. The universal behaviour of the distributions in the number 
of observed jets, just discussed, can be tested in a model independent way, 
which allows a rather strict and detailed control on the different theoretical el- 
ements which enter in the representation of the process. Once the production 
mechanism is tested in detail, possible deviations from the universal behavior 
in particular conditions (heavy nuclei, high energy, central collisions) would 
represent a convincing signal of a different production mechanism taking 
place. 
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Figure captions 



Fig. la. Projection in transverse plane of the scattering centres in a collision 
of a parton with a Pb nucleus. The scale is the fm and the size of the dots 

corresponds to the parton-parton cross section integrated with the cut off in 
momentum transfer A;^ = 5GeV. The momentum fraction xa of the projectile 
parton is xa = 10~^ and the energy in the nucleon-nucleon cm. system is 
7TeV. 

Fig. lb. Same as in fig. la with xa = 10~^. 

Fig. 2a. Rapidity distribution of produced jets in a central Pb — Pb colli- 
sion with 7TeV energy in the nucleon-nucleon cm. system. The histogram 
is the simulated event, the continuous line is the average number of wounded 
partons expected, as from Eq.(7), summing the projectile and target wouded 
partons and expressing the resulting numbers as a function of rapidity. The 
dashed line is the same average number when using the approximate ex- 
pression in Eq.(9), namely after neglecting semi-hard parton rescatterings. 

= 2QGeV. 

Fig. 2b. Same as in fig. 2a with — lOGeV. 
Fig. 2c. Same as in fig. 2a with = 5GeV. 

Fig. 3 a. Distribution in the total number of jets produced in central 
collisions with cutoff k^ = 20GeV. 

Fig. 3b. Same as in fig. 3a with k^ — lOGeV. 

Fig. 3c. Same as in fig. 3 a with k^ = 7GeV. 

Fig. 4a. Relative rates of events with jets aty = 0, k^ = 20GeV in the case 
of a Woods-Saxon distribution in space of the nuclear partonic matter (con- 
tinuous histogram) and for an hard-sphere distribution (dottet histogram). 

Fig. 46. Ratio of the two histograms in fig. 4a. 
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